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Abstract 


Numerical simulations of elementary gravitation and electromagnetic fields are done with an amazingly simple 
algorithm. But this algorithm renders nature correctly. As well, the material world is revealed to be completely 
Riemannian-geometrical without exception. Mathematics is based on the Geometric theorie of fields, which refers 
to Einstein and Rainich. The correctness of the theory is manifested in it that known particles appear as discrete 
solutions of geometric field equations. The results involve new understanding of mathematical principles. 
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1 Introduction 


1 


Physicists commonly believe in four “forces”, if 
not more. Author’s experience of life consists in only 
two, namely gravitation and electromagnetism. Other 
“forces” are tied to specific models introduced to save 
the energy balance of the model. Is a model, that needs 
hidden forces or virtual particles for the energy balance, 
viable at all ? 

The geometric view provides consistent energy bal- 
ance without need to invent additional terms, let alone 
inconsistent models. It is based on General relativity 
[2]. The obstacle consisted in the lack of an algorithm, 
which is suited to naturally solve geometric field equa- 
tions. This gap is filled. 


2 Mathematics 


The complete theory is based on the tensor equa- 
tions [7] 


1 
Rix + 3Ko gik = (F GixFanF” — FiaFk") — (1) 
Fike + Fyn + Fri = 9 ; (2) 
rae =0 ; (3) 


in which g;z are the components of metrics, Riz those 
of the Ricci tensor and Fj, those of the electromagnetic 
field tensor. K, is the constant part of the Riemannian 
curvatures [3], and meaningful for global solutions, e.g. 
[8]. «& is Einstein’s gravitation constant. 

Einstein quoted these equations already in his Four 
lectures [2]. As well, any sources (distributed charges 


1This paper is an enhanced version of a contribution to 
ICMSA-2011, New Delhi [1]. 


and currents) in Eq. (3) would not meet the Bianchi 
identities [3]. The Bianchi identities are mathemati- 
cal expression of force equilibrium respectively energy 
conservation. We shall see that sources are replaced by 
integration constants in the solutions [6]. Mass, spin, 
charge, magnetic momentum are first integration con- 
stants. 

Analytic solutions (different from zero) based on 
integration constants lead commonly to singularities. 
Physicists use to see this as an obstacle. However, nu- 
merical simulations according to Eq. (1),(3),(4) result 
in another picture. [6] 

If we express the field tensor by a vector potential 


Fin = Aik — Ak, : (4) 


Eq. (2) is identically met. Thus, we can go from quan- 
tities with potential character, that are metrics and 
electromagnetic vector potential. 

The geometry of fields is constituted in [6, 7]. Grav- 
itation and electromagnetic fields are parameters of the 
time-like curve described by the test body in the four- 
dimensional space-time. It must be mentioned that 
already Rainich [4, 5] knew the geometry of electro- 
magnetism. This fact seems to be not omnipresent. 
The borderline case of electromagnetic waves (light 
&c.) can be understood from Maxwell’s theory alone 
[9]. 

The geometric equations do not define causality, 
because causality is not a geometric category. Any di- 
alectic philosophy is not sufficient to understand this. 
Smarandache provides an approach to describe this is- 
sue with his “Neutrosophics” [10, 11]. It meets the fact 
that the geometric equations yield only 10 independent 
equations for 14 components giz, A; , considering the 
special role of time. 
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3 Numerical Simulations 


The basical idea consists in it to go from known area 
out of the particle and to successively make toward the 
unknown area. The geometric equations proved in the 
known area, and it is an artificial assumption that these 
be not valid everywhere. It will be demonstrated that 
this artificial assumption is false. 

Basis of the computation are Eq. (1,3,4). For the 
sake of simplicity, we confine it to time independence 
and rotation symmetry. That results, with spherical 
coordinates 

gi=r,2?=0, =, st =jet, 
in 6 independent equations for 8 components A3, Ag, 
G11 9125 922; 933; 934, 944, the rest. vanishes. 
Now, we give 


giz = 0 (with it g'* = 0) (5) 


and 

g = det|giz| = r* sin? 0 . (6) 
That are arbitrary conditions, but they lead to reason- 
able results. 

The integration constants to insert in the initial 
conditions result from development of series. The first 
integration constants are 
Km Km 

= 


= 7 
a An Ar te 
(mass), 
KS KS 
ae 4c 4ic (8) 
(spin), 
2 Q K2 fo? Q 
-_ Lo ° 
ra An Ar (9) 
(charge), and 
1 eg 
Eo2 M K2 €o2 M 
ik Ar At (10) 


(magnetic momentum). 

The dimensionless terms toward the arrow are for com- 
putation, and have positive values. The imaginary unit 
is included. The radius unit (r = 1) is 107-!°m. With 


it, the initial conditions become using T = 5 — J 
= 1628 22 Se Seri ee rh 
gu r 2\r’? ' 10%r?2 : 
C4 1 3 
gn = P{l+(S)*(Feos*T- =}, (12) 
Cc 1 3 
“a = coe Tilt (4) Ge cos” T Te 
CL 1 C3 \2 1 C4\2 . 2 
| T 14 
gua 2 (2)? + (4)? sin (14) 
2 C2 1 c3c4 
gsa = reos T (3 _ rc, ; (15) 
A3 = Reon Ts (16) 
Aye 2 (17) 


The actual computation is done with quantities per- 
formed from physical components, i.e. the unities in 
the equations are eliminated. Because the physical 
components have a magnitude of 1074°. We have to 
insert the values of the integration constants in the 
modified initial conditions, see program in the “robust” 
package (available at author’s current website). Table 1 
shows examples, used as reference values. 


Proton Free electron 
c, | 2.48 x 10-89 | 1.30 x 10747 
c2 | 2.60 x 1074° | 2.60 x 10749 
cz | 1.95 x 1077! | 1.95 x 107?! 
C4 5.7 x 107-72 3.7 x 107-19 

Deuteron Helium nucl. 
c, | 4.96 x 107%9 9.9 x 107°9 
co | 5.2 x 10-40 0 
cz | 1.95 x 107?! 3.9 x 10721 
c4 | 1.76 x 10722 0 

Table 1 


As well, the higher momenta are missing for lack of 
knowledge, what affects the correctness of the compu- 
tation. We will insert known values and values deviat- 
ing from them, and compare the results. 

The algorithm consists in performing difference 
equations, in which the newly to quantifying quantity 
is on the left-hand side, and all previously quantified 
quantities are on the right-hand side. If we calculate 
spherical shells from outside to inside, the new quantity 
is fm+2,n in following difference equations (f stands for 
any potential quantity) 
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Of 7 fate _ Stan 

or Tmln ~ 2dr : (18) 
O° f = fm+2.n ss 2fmin -F fim—2,n (19) 
Or? |r Tn (2 dr)? 

Of = fmm41 = eel 

OT lrm;Tn 2dr (20) 
oF fmn+1 = 2fmin ak Tw 
Wn, = ar” oe) 


This quantity is linearly dissolvable in Eq. (19), and 
the whole adherent tensor equation is non-linear. De- 
tailed formulae are to see in the Pascal codes. (The 
current code is in the “robust” package.) 

When the program runs, the values of the several 
components are successively quantified in one spherical 
shell after the other (from previous four outer). We get 
at first the values as expected from the initial condi- 
tions. That runs so a while. Suddenly, the values grow 
out all limits. The computation is broken when geo- 
metric limits are reached. — The tendency is following: 
If we halve the differences, we get noway double num- 
ber of steps (until appearance of this effect) but clearly 
less. That would mean: There are no solutions at all. 
But there are exceptions for certain values of the inte- 
gration constants ! — That has to do with chaos, and 
it is the job of mathematicians to class this. 

In order to see these exceptions, lots of tests were 
done with parameters more and less deviating from the 
references (Table 1). Output is mentioned number of 
steps. The differences of the step numbers are small but 
correlate with the known quantities of particles. As 
well, the correlations became highly significant when 
the raster distances were the same tangentially as well 
as radially (dr = r d@ ) just at the conjectural particle 
radius. 

Fig. 1 and 2 show visualized results from these 
tests. The number of steps above a “threshold” is de- 
picted with a more or less fat “point”. 


4 Conclusion 


We can understand nature completely Riemannian- 
geometrically as far as we own to the fact that nature 
follows neither analytic solutions nor models invented 
to explain special effects. A natural solution of par- 
tial differential equations goes always from finite dif- 
ferences. That leads to a new understanding of partial 
differential equations, and the Leibniz calculus in gen- 
eral: The transition to differences towards zero may be 
done first at the end of all calculations. But Konrad 
Zuse [12] asked the question if this transition is pos- 
sible at all. However, a discretized space-time suffers 


from the problem of preferred coordinates. These were 
not verified in nature. 
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y-axis: charge 1.0E-21 to 9.1E-21 Helium nucleus, small circle = reference, 82A2 tests 


Mass of: {proton fdeuteron tHelium nucleus x-axis: mass 2.0E-39 to 2.225E-38 


Figure 1: Tests with parameters around the Helium nucleus 


y-axis: charge 1.0E-21 to 3.0E-21 Electron, small circle = reference, 41A3 tests 


x-axis: magnetic momentum 1.0E-19 to 5.0E-19 , perspective: spin 1.7E-40 to 3.7E-40 


Figure 2: Tests with parameters around the electron 
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